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Forced turbulence combined with the effect of rotation and shear flow is studied. In a previous 
paper [Leprovost and Kim, PRE in press (2008)], we considered the case where the shear and the 
rotation are perpendicular. Here, we consider the complementary case of parallel rotation and shear, 
elucidating how rotation and flow shear influence the generation of shear flow (e.g. the direction of 
energy cascade), turbulence level, transport of particles and momentum. We show that turbulence 
amplitude and transport are always quenched due to strong shear (£ = uky/A <C I, where A is 
qq . the shearing rate, v is the molecular viscosity and k y is a characteristic wave-number of small-scale 

turbulence), with stronger reduction in the direction of the shear than those in the perpendicular 
directions. In contrast with the case where rotation and shear are perpendicular, we found that 
£NJ ■ rotation affects turbulence amplitude only for very rapid rotation (fi 3> A) where it reduces slightly 

! > ' the anisotropy due to shear flow. Also, concerning the transport properties of turbulence, we find 

O \ that rotation affects only the transport of particle and only for rapid rotation, leading to an almost 

isotropic transport (whereas, in the case of perpendicular rotation and shear, rotation favors isotropic 
transport even for slow rotation). Furthermore, the interaction between the shear and the rotation 
is shown to give rise to non-diffusive flux of angular momentum (A-effect), even in the absence of 
external sources of anisotropy, which can provide a mechanism for the creation of shearing structures 
(— i ■ in astrophysical and geophysical systems. 

>v 

^ , PACS numbers: 47.27. Jv,47.27.T-,97.10.Kc 

i : 

I. INTRODUCTION 

Large-scale shear flows are often observed in rotating astrophysical and geophysical systems. Shear and rotation 
have a huge impact on the properties of this system, such as energy transfer or mixing. From a physical point of view, 
the main effect of shear flow is to reduce turbulence level as well as turbulent transport compared to their values 
without shear. This is basically because shear advects turbulent eddies differentially, elongating and distorting their 
shapes, thereby rapidly generating small scales which are ultimately disrupted by molecular dissipation on small scales 
(see Fig. [T|). That is, flow shear facilitates the cascade of various quantities such as energy or mean square scalar 
density to small scales (i.e. direct cascade) in the system, enhancing their dissipation rate. As a result, turbulence 
level as well as turbulent transport of these quantities can be significantly reduced compared to the case without shear. 
' Another important consequence of shearing is to induce anisotropic transport and turbulent level since flow shear 
y—( , directly influences the component parallel to itself (i.e. x component in Fig. 1) via elongation while only indirectly 
£N| ' the other two components (i.e. y and z components in Fig. 1) through enhanced dissipation. Rotation can also 
• reduce transport in the limit of rapid rotation (similarly to flow shear), but through a physical mechanism that is 
different from that of shear, namely by phase mixing of inertial waves While both rotation and (stable) shear flow 
tend to regulate turbulence, there are important differences in their effects, which should be emphasized. Rotation, 
by exciting inertial waves, tends to reduce turbulence transport more heavily than turbulence amplitude while shear 
flows reduce both of them to a similar degree. That is, rotation (or waves) quenches the cross-phase (normalized flux) 
• i-H . more than shear flow does 0, Q . 

Rapid distortion theory (RDT) 0, [j| was used to study the linear response of turbulence to a mean flow with 
spatially uniform gradients. The linear treatment of fluctuations by incorporating strong flow shear was also used in 
the astrophysical context by @ by using shearing coordinates. The generation of large-scale shear flows (the so-called 
zonal flows) through a similar nonlocal interaction has been intensely studied in the magnetically confined plasmas, 
where turbulence quenching by shear flow is believed to be one of the most promising mechanisms for improving 
plasma confinement 0,0- In decaying sheared turbulence, Q have shown a surprisingly good agreement between the 
RDT predictions and numerical simulations. Subsequently, theoreticalpredictions (using a quasi- linear theory) for the 
transport of passive scalar fields in 2D hydrodynamic turbulence by Q and [jfj| have been beautifully confirmed by 
recent numerical simulations [Til ] . In particular, they have shown that turbulent transport of particles can be severely 
quenched inversely proportional to flow shear A while turbulence level is reduced as A~ 5 ^ 3 . It is important to note that 
this nonlocal interaction leading to inverse cascade can be successfully captured by inhomogeneous RDT theory which 
permits the feedback of the nonlinear local interaction between small scales onto the large scales via Reynolds stress 
(constituting the other part of quasi-linear analysis) while neglecting nonlinear local interaction between small scales 
for fluctuations compared to nonlocal interactions. As must be obvious by comparing the Coriolis force with nonlinear 
advection terms, the RDT works well for sufficiently strong rotation (small Rossby number) even in the absence of 
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FIG. 1: Sketch of the effect of shear on a turbulent eddy. 



shear flow. For instance, the agreement of the RDT prediction with numerical results has been shown by various 
previous authors including [r3 |. but mostly in decaying turbulence. However, in this case, the RDT cannot accurately 
capture the turbulence structure in the plane perpendicular to rotation axis where nonlinear local interactions between 
inertial waves seem important (see, e.g. [l3]). The validity and weakness of the RDT together with comparison with 
various numerical simulation (without an external forcing) with/ without shear flows and stratification can be found 
in excellent review by ,14] and Cambon and [15l ]. to which readers are referred for more details. 

In comparison, far much less is understood in the case of forced turbulence. In particular, the main interest 
in forced turbulence is a long-term behavior where the dissipation, enhanced by shear distortion, is balanced by 
energy input, thereby playing a crucial role in leading to a steady equilibrium state. The computational study of 
this long time behavior is however not only expensive but also difficult because of the limit on numerical accuracy, 
as noted by (l6| . Therefore, analytical theory by capturing shearing effect (such as quasi-linear theory with time- 
dependent wavenumber) would be extremely useful in obtaining physical insights into the problem as well as guiding 
future com put ational investigations. We note that the previous works by Kichatinov and Rudiger and collaborators 
[13, d, fl9l. l2fj|. l2ll ] using quasi-linear theory are valid only in the limit of weak shear. Forced sheared turbulence was 
proposed for the first time by 22[ in the context of two-dimensional near-wall turbulence to explain the logarithmic 
dependence of the large scale velocity on the distance to the wall. In that case, the external forcing is provided by 
a continuous supply of vorticity from intermittent coherent burst of vorticity coming from the viscous layer. This 
work was later generalized to three dimensions [23l |24| with the same conclusions. In the astrophysical context, 
[25| has shown that in 3D forced HD turbulence, strong flow shear can quench turbulence level and transport of 
particles with strong anisotropy (much weaker along the flow shear which is directly affected by shearing) and has 
emphasized the difference in turbulence level and transport, which is often used interchangeably in literature. A 
similar weak anisotropic transport was shown for momentum transport by [2o| in forced 3D HD turbulence. Further 
investigations have been performed on turbulent transport in forced turbulence by incorporating the interaction of 
sheared turbulence with different types of waves that can be excited due to magnetic fields 0, H3, ,28j] , stratification 
[29l | or both magnetic fields and stratification (30} . 

The combined influence of shear flow and rotation in forced turbulence has been considered by Leprovost and Kim 
in paper I [3l| when the rotation and the shear are perpendicular to each other. We found that flow shear always 
leads to weak turbulence with an effectively stronger turbulence in the plane perpendicular to shear than in the shear 
direction, regardless of rotation rate. The anisotropy in turbulence amplitude is however weaker in the rapid rotation 
limit (fi 3> A where fi and A are the rotation and shearing rate) than that in the weak rotation limit (il <C A) 
since rotation favors almost-isotropic turbulence. Compared to turbulence amplitude, particle transport is found to 
crucially depend on whether rotation is stronger or weaker than flow shear. When rotation is stronger than flow shear, 
the transport is inhibited by inertial waves, being quenched inversely proportional to the rotation rate while in the 
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FIG. 2: Sketch of the configuration in the parallel case 



opposite case, it is reduced by shearing. Furthermore, the anisotropy is found to be very weak in the strong rotation 
limit (by a factor of 2) while significant in the strong shear limit. The turbulent viscosity is found to be negative 
with inverse cascade of energy as long as rotation is sufficiently strong compared to flow shear while positive in the 
opposite limit of weak rotation . Even if the eddy viscosity is negative for strong rotation, flow shear, which transfers 
energy to small scales, has an interesting effect by slowing down the rate of inverse cascade with the value of negative 
eddy viscosity decreasing as |z/y| oc A~ 2 for strong shear. Furthermore, the interaction between the shear and the 
rotation is shown to give rise to a novel non-diffusive flux of angular momentum known as the anisotropic kinetic 
a-effect (AKA) [32[ or as the A-effect in the astrophysical community. The appearance of non-diffusive term in the 
transport of angular momentum prevents a solid body rotation from being a solution of the Reynolds equation [33L [HH , 
and thus act as a source for the generation of large-scale shear flows. For instance, this effect has been advocated 
as a robust mechanism to explain the differential rotation in the solar convective zone. Starting from Navier-Stokes 
equation, it is possible to show that these fluxes arise when there is a cause of anisotropy in the system, either due 
to an anisotropic background turbulence (see [20| and references therein) or else due to inhomogeneities such as an 
underlying stratification. In [3l|, we found that a A-effect appears in sheared-rotating turbulence even in the absence 
of external sources of anisotropy. This is because the shear induces an anisotropic turbulence which combined to the 
rotation gives rise to non-diffusive fluxes. 

In this paper, we consider the complementary case when rotation and shear are parallel to each other. By assuming 
either sufficiently strong shear or rotation rate, we employ a quasi-linear analysis to compute turbulence level, eddy 
viscosity, and particle transport for temporally short-correlated, homogeneous forcing. As the computation of these 
quantities involve too complex integrals to be analytically tractable, they are analytically computed by assuming 
an ordering in time scales. In our problem, there are three important (inverse) time-scales: the shearing rate A, 
the rotation rate and the diffusion rate T> = vk 2 , where v is the (molecular) viscosity of the fluid and k~ l is a 
characteristic small scale of the system. We first distinguish the two cases of strong rotation (f2 ^>> A) and weak 
rotation (fl A). The first regime of strong rotation will be studied in the strong shear (A ^> T>) and weak shear 
{A <C T>) regime. On the other hand, the second regime of weak rotation will be considered only in the strong shear 
{A T>) case, as the effects of both shear and rotation disappear in the opposite limit (icP). 



II. MODEL 

We consider an incompressible fluid in a rotating frame with average rotation rate f2, which are governed by 

ftu + u-Vu = -VP + zA7 2 u + F - 2fi x u , (1) 
V ■ u = . 

Following [251 ]. we study the effect of a large-scale shear U = Uo(x)j on the transport properties of turbulence by 
writing the velocity as a sum of a shear (chosen in the x-direction) and fluctuations: u = Uq + v = Uo(x)j + v = 
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—xAj + v. Without loss of generality, we assume A > 0. In the following, we consider the configuration of Figure 
[5] where the shear and rotation (in the x direction) are parallel and simplify notation by using O = 2J~2. Then, the 
Coriolis force is simply f2[— v z j + v y k], where i, j and k are the unit vectors associated with the Cartesian coordinates. 
Note that our x — y coordinates are not conventional in that our x and y directions correspond to y and x in previous 
works (see [l6j for instance). Therefore, the shearing, the stream- wise and the span- wise direction correspond to the 
x, y and z direction, respectively. Note that the particular configuration considered here, the stationary equation for 
the large-scale velocity reduces to: 

-d x P + F x = 0, (2) 

-d z P + F z - nu Q (x) = 0. 

For F = 0, there is no equilibrium for the large-scale flow contrary to the perpendicular case where the Coriolis force 
can be balanced by pressure. Thus, we assume a large-scale forcing F to maintain an equilibrium for the large-scale 
flow [e.g. [35| . In real physical situation, this lar ge-s cale forcing can be provided by thermal wind associated with 
(latitudinal) temperature gradient on large scales [Til l3o| . Density/temperature fluctuations in this case can then be 
included as a part of the small-scale external forcing f in our formalism. Note that a similar approach was adopted by 
[3(| in computing the transport of momentum and particle in stars, including the polar region. A consistent treatment 
of fluctuating temperature requires the extension of the present work to stratified rotating sheared turbulence and 
will be performed in future publications. 

To calculate turbulence amplitude (or kinetic energy) and turbulent transport, we need to solve the equation for 
the fluctuating velocity field. To this end, we employ a quasi-linear theory [371 ] where the nonlinear local interactions 
between small scales are neglected compared to nonlocal interactions between large and small scales and obtain: 

d t v + U ■ Vv + v • VU = -Vp + ;A7 2 v + f - fl x v , (3) 
V • v = , 

where p and f are respectively the small-scale components of the pressure and forcing. As noted in the introduction, 
this approximation, also known as the RDT :5j, is justified in the case of strong shear as the latter induces a weak 
turbulence, leading to weak interaction between small scales which is negligible compared to the (non-local) interaction 
between the shear and small scales. This has in fact been confirmed by direct numerical simulations, proving the 
validity of the predictions of quasi-linear theory with a constant-rate shear both in the non-rotating 9] and rotating 
unforced [l6[ turbulence and also for forced turbulence [llj. Further, note that the quasi-linear analysis is also valid 
in the limit of rapid rotation [38| . 

To solve Eq. ([3]) , we introduce a Fourier transform with a wave number in the x direction evolving in time in order 
to incorporate non-perturbatively the effect of the advection by the mean shear flow 0, [^, [25| : 

v (x,£) = _I_ J d 3 k e l t k ^ x+k «y +k * z ^(k,t) , (4) 

where k x (t) = k x (0) + k y At. From Eqs. ((3j) and j4|), we obtain the following set of equations for the fluctuating 
velocity: 

Ad T v x = -ik y rp + f x , (5) 

A8 T Vy - AV X = ~ikyf> + fy + flV Z , 

Ad T v z = -ik z p + f z - Qv y , 
= tv x +v y + (3v z . 

Here, the new variables v = v exp[^(fc| f t + k^/3k y A)) and similarly for f and p have been used to absorb the diffusive 
term, and the time variable has been changed to r = k x (t)/k y . In the remainder of the paper, we solve Eq. ([5]) for 
the fluctuating velocity (with a vanishing velocity as initial condition) . We then use these results and the correlation 
of the forcing (defined in mi CP to compute the turbulence intensity and transport (defined in mi B[) . 

A. Transport of angular momentum 

As the large-scale velocity is in the y direction, we are mostly interested in the transport in that direction. The 
large-scale equation for the y component of velocity Uo is given by Eq. |T]) with a supplementary term V • R where 
R is the Reynolds stress given by: 

R = (vv v ) . (6) 
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To understand the effect of R on the transport of angular momentum, one can formally Taylor expand it with 
respect to the gradient of the large-scale flow: 

Ri = AiU - v T d x U 5n + ■ ■ ■ = A,U + v T A8n + ... . (7) 

Here, Aj and vt are the two turbulent transport coefficients from non-diffusive and diffusive momentum flux, re- 
spectively. Note that the first term in the expansion is due to the small-scale driving and the Coriolis force in Eq. 
§\§ which break the Galilean invariance [39 ]. First, vt is the turbulent (eddy) viscosity, which simply changes the 
viscosity from the molecular value v to the effective value v + vt- Note that the sign of eddy viscosity represents the 
direction of energy cascade, with positive (negative) value for direct (inverse) cascade. Second, the first term involving 
Ai in equation ([7]) is proportional to the rotation rate rather than the velocity gradient. This means that it docs not 
vanish for a constant velocity field and thus permits the creation of gradient in the large-scale velocity field. This 
term bears some similarity with the a effect in dynamo theory [5(| El[ and has been known as the A-effect [13, 
or anisotropic kinetic alpha (AKA)-effect [12]. Similarly to the a effect, this effect exists only if the small-scale flow 
lacks parity invariance (going from right-handed to left handed coordinates). However, in contrast to the a effect, the 
A effect requires anisotropy for its existence [13, HH . 



B. Particle (or heat) transport 

To study the influence of rotation and shear on the particle and heat transport, we have to supplement equation ([T|) 
with an advection-diffusion equation for these quantities. We here focus on the transport of particles since a similar 
result also holds for the heat transport. The density of particles N(x,t) is governed by the following equation: 

d t N + U • VN = DV 2 N , (8) 

where D is the molecular diffusivity of particle. Note that, in the case of heat equation, D should be replaced by 
the molecular heat conductivity x- Writing the density as the sum of a large-scale component No and small-scale 
fluctuations n (N — No + n), we can express the evolution of the transport of chemicals on large scales by: 

d t N + U • VJV„ = (D5 i:i + Dij> )d i d j N , (9) 

where the turbulent diffusivity is defined as (w,n) = —Dj,djNo. The turbulent diffusivity is computed in the following 
to study the effect of rotation and flow shear on turbulent transport of chemicals which can be highly anisotropic. 
Note that the transport of a passive scalar quantity (contrary to the angular momentum which is a vector quantity) 
has to be diffusive due to the fact that it is solely advected by the flow [42]. 

For simplicity, we assume a unit Prandtl number D = v and apply the transformation introduced in equation ((3]) 
to the density fluctuation n to obtain the following equation: 

d T n = HW fl . (10) 
A 

Equation (|10p simply shows that the fluctuating density of particles can be obtained by integrating the fluctuating 
velocity in time. 



C. External forcing 



As mentioned in introduction, we consider a turbulence driven by an external forcing f . To calculate the turbulence 
amplitude and transport defined in §11 Al and flll Bl (which involve quadratic functions of velocity and/or density), we 
prescribe this forcing to be short correlated in time (modeled by a (^-function) and homogeneous in space with power 
spectrum tpij in the Fourier space. Specifically, we assume: 

(fi(^l,ti)fj(^2,t 2 )) =r / (2^) 3 «5(k 1 +k 2 )«5(t 1 -i 2 )^(k 2 ) , (11) 

for i and j = 1, 2 or 3. The angular brackets stand for an average over realizations of the forcing, and r/ is the (short) 
correlation time of the forcing. Note that the 5 correlation is valid as long as the correction time 17 is the shortest 
time-scale in the system [i.e. 77 <C £l~ l ,A~ l , \/{vk 2 )]. 



6 



For most results that will be derived later, we assume an incompressible and isotropic forcing where the spectrum 
of the forcing is given by: 

^ ij (\ S )=F(k)(8 ij -k i k j /k 2 ) . (12) 
It is easy to check that in the absence of rotation and shear, this forcing leads to an isotropic turbulence with intensity: 

where the subscript stands for a turbulence without shear and rotation. 



III. ANALYTICAL RESULTS 



To investigate turbulence property in the parallel case, we rearrange ([5]), thereby eliminating the pressure terms, 
obtaining the following equations for fluctuating velocity: 
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dr{{l + T 2 }v x ) 
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d T Vz = 
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(14) 
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Here, (3 = k z /k y , 7 = 
equation in equation 
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k 2 H jk 2 and Q 



1 A. To calculate the turbulence amplitude and transport, the first 



1 

14]) is to be solved with initial conditions: v x (tq) — and d T v x \ T=To = /ii(tq)/(7 + t 2 )A. 
Unfortunately, we were unable to find an exact solution of equation (|14|) in the general case. To gain a physical 
insight into the role of inertial waves and flow shear in turbulent transport, we consider the two limits of strong 
rotation where the effect of waves dominates shearing (51 3> A) and weak rotation where shearing dominates the 
effects of waves (51 <C A) in section IIII Al and IIII Bl respectively. Approximate solutions can be derived in these 
two regimes which can then be used for deriving analytic form of correlation functions for turbulence intensity and 
transport. However, the WKB approximation used to study the rapid rotation limit being ill-defined for some value 
of the parameter, we consider an exactly solvable case where all quantities can be calculated in section UlI CI 



A. Rapid rotation limit: Q 2> A (and a > 0) 



For |0| 3> 1, we seek for a WKB solution of the first equation in equation (|T4"1) . However, since this approximation 
breaks for r = 0, we assume that the initial value is positive (tq = a > 0) to make our solution meaningful. In mil CI 
we study an exactly solvable case and show that the solution can be altered by negative initial value (tq — a < 0). 
Assuming tq = a > 0, we obtain the following solutions for the three components of the velocity for |51/,4| ^> 1: 
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Due to the similarity between equation (| 1 5|) and the one obtained in the perpendicular case [3l| , we give only results 
and refer the reader to our previous paper for details of the derivation. 



1. Turbulence intensity 

In the weak shear limit (£ 3> 1) where the shear is negligible, we obtain the following result for the turbulent 
intensity: 



VI 



3(2tt) 



Here, (vq) is the turbulence amplitude in the absence of rotation and shear [see Eq. ([T3]) ]. These results thus show 
that, in the large rotation limit, the turbulence intensity is isotropic and equals to the one without rotation [see Eq. 
(flBf ] for sufficiently weak shear with £ 3> 1. Furthermore, in this limit of sufficiently weak shear where (ft,T>) ^s> A, 
turbulence intensity is independent of rotation since waves do not necessarily quench turbulence level. A similar result 
was also obtained in MHD turbulence and stratified turbulence where magnetic fields and gravity waves mainly affect 
transport without much effect on turbulence level [H, H^, [3(j • We shall show below that a strong anisotropy can be 
induced when shearing effect is not negligible (£ -C 1) even in the rapid rotation limit (SI ^> .4). 

In order to understand the effect of flow shear, we now consider the strong shear limit (( < 1). In this limit, we 
obtain: 

(«*) = j^fA J Sk ^ + a2p W cc £(« 2 ) , (18) 



yl " z/ (2tt)M 



d 3 k y^Ta^F(k)—^ cx £|ln£|(t, 2 ) , 



to leading order in £ <C 1. Note that in the calculation of (t> 2 ), we neglected the component proportional to a = k x /k y 
as it is odd in both k x and k y and thus vanishes after integration over the angular variables for an isotropic forcing. 
The last terms in Eq. ((T5J), expressed in terms of the turbulence amplitude in the absence of rotation and shear (v^) 
[see Eq. (|13p]. explicitly show the dependence of turbulence level on rotation and shear. That is, all the components 
of turbulence intensity are reduced for strong shear ( < 1. Further, the x component along shear is reduced as 
£ oc A^ 1 while the other two components as £| ln£|, with an effectively weaker turbulence in the shear direction than 
in the perpendicular one, by a factor of | m£|. This shows that shear flow can induce anisotropic turbulence (unlike 
rotation) even when the forcing is isotropic. This result is similar to that obtained in the simulation of a Couette 
flow at high rotation rate [H[ where the velocity fluctuations perpendicular to the wall exceed that in the stream- wise 
direction. Nevertheless, Eq. (| 18() shows that a strong rapid rotation yet insures an isotropy in velocity fluctuations in 
y — z directions ((i; 2 ) = (« 2 )). 



2. Transport of angular momentum 

First, in the large rotation and weak shear limit, the Reynolds stress becomes purely diffusive (with no A-cffcct) 
with the turbulent viscosity: 

7TT f f + °° „ F(k) , , 

VT — ^_ / dk— — . (19) 

This result shows that the turbulent viscosity is positive and proportional to SI -1 for large fi. 

In comparison, in the strong shear limit (£ <C 1), we obtain the turbulent viscosity in the strong shear limit as: 

-=^=-<^W af «- < 2 »> 

Eq. ([2H)l shows that the turbulent viscosity is negative (as F(k) > 0) in the strong shear limit, in sharp contrast 
to the weak shear limit where vt > [see Eq. (|19l) ]. Furthermore, the magnitude of Vt is reduced by the shear 
(cx A~ 2 ) and is independent of rotation, which should also be compared with the weak shear limit [see Eq. (TH?|) where 
Vt oc SI -1 ]. Therefore, the turbulent viscosity changes from positive (for weak shear) to negative (for large shear) as 
the ratio of shear to dissipation increases. This result can be understood if we assume that, as in most rapidly rotating 
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fluid, the inverse cascade is associated with the conservation of a potential vorticity [44]. In the presence of strong 
shear (compared to dissipation), the potential vorticity is strictly conserved giving rise to an inverse cascade (negative 
viscosity). When the dissipation increases, the potential vorticity is less and less conserved and thus the inverse 
cascade is quenched. Our results show that there is a transition from inverse to direct cascade as the dissipation is 
increased. A similar behavior is also found in two-dimensional hydrodynamics (HD) where an inverse cascade can be 
shown to be present only for sufficient weak dissipation (27j |. 

It is important to note that the negative viscosity vt < obtained here for strong rotation/strong shear (fi 3> 
A vky) signifies the amplification of shear flow as the effect of rotation favoring inverse cascade dominates shearing 
(generating small scales). However, the magnitude of vt is reduced by shear as \vt\ c< A~ 2 since flow shear inhibits 
the inverse cascade. This can be viewed as 'self-regulation' - that is, self-amplification of shear flow is slowed down 
as the latter becomes stronger. 



3. Transport of particles 
In the rapid rotation limit (|f2|/„4 3> 1), turbulent particle diffusivities can be obtained as: 

d t ~ r — ^ > ( 2i ) 

yy _ nzz _ T f f jl „, 



D™ = D Z T Z ~ — J — / dk ~ -D 

T T 16tt|0| J q v 2 ^ 

Note that in that case, the result is not sensitive to the value of the parameter £ and thus we do not distinguish 
between the weak and large shear limits. Eq. (121"]) shows that D^ x , Dip 1 and Dff are all reduced as f2 _1 (with no 
effect of the shear) for large and also that there is only a slight anisotropy in the transport of scalar: the transport 
in the direction of the rotation is twice larger than the one in the perpendicular direction [2l| . Interestingly, this 
anisotropy in the transport of particles is not present in turbulence intensity [see Eq. JT7J]. This is because waves 
can affect the phase between density fluctuation and velocity, not necessarily altering their amplitude. However, it 
is important to note that this anisotropy is only a factor of 2, much weaker than that in sheared turbulence without 
rotation 12511. 



B. Weak rotation limit: Q <C A 



In the weak rotation limit, we expand all the quantities in powers of flo = \Q/A\ as: 

X(T)=X ( T )+n X 1 (T) + ... , (22) 

in the weak rotation limit (Q ^ ^4) and calculate the turbulence intensity and transport up to first order in f^o- 
For the sake of brevity, we here just provide the final results of the calculation. Note that in this limit, we are only 
interested in strong shear case (£ 1) since in the opposite limit where vk 2 ^S> A ^S> ft, the effects of both shear and 
rotation simply disappear to leading order. 



1. Turbulence intensity 



In the strong shear limit (£ <§; 1), we obtain the turbulence intensity up to order as follows: 
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d 3 kF(k) 



J . .2 ( 71 
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-T(a) +1 
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V2f 



1/3 



1/3 



r(i/3). 



(23) 



Here, T is the Gamma function. Note that the first correction (proportional to 17) vanishes and consequently, the 
turbulence amplitude is the same as in the case of shear without rotation [2 51 ] , with stronger turbulence in the direction 
perpendicular to the shear than in the parallel one. 
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2. Transport of angular momentum 



In the strong shear limit (£ <C 1), momentum flux in the azimuthal direction can be shown to be purely diffusive 
and given by: 



{V X Vy) 



(2tt)M 



d 3 k { 1 + a 2 )F{k) 



1 



2( 7 + a 2 ) 



■a 2 



2^7 



-T(o) 



(24) 



This recovers the eddy viscosity of sheared turbulence without rotation [25|, showing that its value decreases as 
oc A~ 2 for strong shear. This result agrees with previous studies of non-rotating sheared turbulence [23| which found 
a Reynolds stress inversely proportional to the shear, leading to a log dependence on the distance to the wall for the 
large-scale shear flow. 

Alternatively, the component of the Reynolds stress involving the velocity component v z does not vanish and is odd 
in Q. Thus, the A-effect appears here in the z-component of the Reynolds stress A z (recall that in the perpendicular 
case, the A-cffcct was present only in (v x v y )), in the form: 



A, 



(2tt)M 2 



d 3 kF(k) 



£(2/3) 



|) 2/3 (3^-l){( 7 + fl ^(^ 



T{a) 



(25) 



Equation ([25]) shows that the sign of A z is indefinite (as both signs appear in the prefactor (3/3 2 — 1). However, as for 
equation (|24p . in the case of an isotropic forcing, the term proportional to 1 dominates, making A z negative. This 
A effect appears due to the anisotropy between the stream-wise and the span-wise components of the velocity, due to 
the shear (alone). 



3. Transport of particles 
Up to order f2, we find the turbulent diffusivity of particles as: 



D X T X 



{2n) 3 A 2 

T f 
(2ir) 3 A 2 



d 3 k~f(j + a 2 )F(k) 



( n 



V2V7 



d 3 k F(k) 



Fh + a 2 ) 



V2V7 



T(a] 



-T(a) +1 



(26) 



^(l) 2,3r(2/3) 



Here again, the first correction due to rotation vanishes and one recover the result of turbulence in presence of shear 
alone: D^ x oc A~ 2 and D^? oc A~ 4 ^ 3 , with effectively faster transport in span- wise direction compared to shear 
direction. 

C. Symmetric perturbation (f3 = 0) 

In this section, we consider a symmetric perturbation with k z = by assuming a forcing that is symmetric in the 
span-wise direction with no dependence on z. Note that even though k z — 0, v z and v y are closely linked through 
rotation fii. Details of the derivation are given in Appendix [X] The interest of this case is that we obtain solutions 
for arbitrary values of f2o and so we can look at the f2o 3> 1 limit without the a > limit. Consequently, we will show 
here only the results in the large rotation limit. 



1. Turbulence amplitude 

In the large rotation limit: f^o 3> 1, we obtain the following leading order contribution of the turbulent amplitude: 

/..2\ _ T f 



(2ir)M 



d 3 kF(k)y/l + a 2 , 



(27) 
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Thus, the turbulence amplitude is larger in the y — z plane than the one in the shear direction by a logarithmic 
factor. Moreover, equation (|27[) shows that the turbulence amplitude does not depend on the rotation rate in the 
large rotation limit but is quenched by shear only. In particular, (Vy) = (v 2 ). These results are the same as in the 
case where the shear and the rotation are perpendicular [3l[ and thus agree with the WKB solution of Sec. IIII Al 



2. Turbulent transport of momentum 



In the large rotation limit (f^o 3> 1) ; we obtain the following turbulent viscosity: 



VT = -j2^ j SkF ^- ( 28 ) 

Equation (|28|) shows that the turbulent viscosity does not depend on rotation in the large rotation limit and is 
obviously negative. Note that this result is the same as in the perpendicular case [see equation (|20[) ] and, thus again, 
agrees with the WKB solution found previously. 



3. Particles transport 



In the limits of strong shear (£ <C 1) and large rotation (f2o ^> 1), the transport of particles is given by: 

d ? ~ mir^ dk+ wiL d3kVTT ^ Fik) - 

The transport of particles in equation involves two contributions, both of which scale as f2 _1 for rapid rotation. 
The first contribution comes from the integration by parts and has to be kept only because ZJq can vanish for a — 
while the second comes from the stationary point in the integration (see appendix IB 31 for details). Note that the ratio 
of the second term to the first one is equal to vk 2 / A ~ £. Consequently, in the strong shear limit (£ <C 1), the first 
term dominates. Thus, the transport of particles is the same as the one found with the WKB analysis (see S jlll Al) . 

To summarize, in this section, we solved equation (fl4|) exactly for j3 = and compared the results with the WKB 
analysis performed in mil Al f which is valid only for a > 0). The results being the same, the conclusions reached from 
WKB analysis remain valid even if a < 0. 



IV. CONCLUSION 



In this paper, we have performed a thorough investigation of the combined effects of shear and rotation on the 
structure of turbulence, by using a quasi-linear theory. We assumed an external forcing in the Navier-Stokes equation 
which leads to an equilibrium situation where the dissipation (whose effect is enhanced by the shear) is balanced by 
the injection of energy due to forcing. It is useful to recall that there are three (inverse) time-scales in the problem: 
the shearing rate A, the rotation rate and the diffusion rate T> — vk 2 where v is the (molecular) viscosity of the 
fluid and l/k y is a characteristic small-scale of the forcing. The first regime of strong rotation (fi 3> ^4) has been 
studied in the strong shear [A ^> T>) and weak shear (.4 <C T>) limits. However, the second regime of weak rotation 
has been considered only in the strong shear (^4 ^> T>) case, as the effects of both shear and rotation disappear in the 
opposite case. 

In the large rotation limit (|fi| ^> A), we found that the results coincides with that obtained in the case where the 
rotation and the shear are perpendicular. Specifically, we obtained the following results: 

• The turbulent intensity is reduced only by a strong shear (i.e. in the case of strong rotation and strong shear) 
and in an anisotropic way. 

• As the dissipation decreases (compared to the shear), there is a crossover from a positive to a negative viscosity. 

• The transport of particle is reduced by rotation, with a slight anisotropy of a factor 2, largely unaffected by 
shear. 
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Perpendicular [311 


Parallel 




n > A 


n < -4 


n > A 


< .4 




A' 1 


^r 1 [l + cn] 


A' 1 






A-^k^] 


_4-2/3 [i + C -n|lnC|] 


•A _1 |ln£| 


^-2/3 




-A' 2 


yU 2 


-A' 2 


^ 2 







^- 2 |lne| 








A z 











A- 4/3 




fr 1 


.4" 2 [l + CQ|lnf|] 




A- 2 




n- 1 


A -4/s [i + C -n|ln^|] 


fi" 1 


A- A ' 3 



TABLE I: Summary of our results obtained for the perpendicular [3l[ and parallel cases in the strong shear limit (£ = vky/A -C 
1) and for an isotropic forcing. In the perpendicular case, the rotation is in the z direction whereas it is in the x direction in 
the parallel case. In both cases, the shear is in the x direction. The C symbol stands for an additional constant of order 1. 

In the opposite weak rotation limit (|0| -C A), we found that the main reduction is due to the shear with an 
anisotropic turbulence with preferred motion and transport in the plane perpendicular to the shear. Contrary to the 
perpendicular case, we found here that the turbulence intensity and the particle transport are not affected by rotation. 
Furthermore, we found non-diffusive flux for momentum transport (the so-called A-effect) which transfers energy from 
the fluctuating velocity field to the large-scale flow. This effect can appear even for isotropic forcing due to the fact 
that the shear induces an anisotropic turbulence. In this paper, the lambda effect appear on the z component which 
has to be contrasted with the perpendicular case [31] where the non-diffusive fluxes appeared on the x-component. 

Table U summarizes the findings of this paper together with these of paper II [3l[ by highlighting the quenching 
of these quantities due to large shearing rate A and the rotation rate £1 (or their ratio, fl = ft /A). We choose to 
show the result in the strong shear limit (£ = vky/A <C 1) as it is the proper limit to capture the effect of the shear. 
Furthermore, from the physical point of view, it is the meaningful limit in a vast number of systems (for example the 
Sun). 

These results can have significant implications for astrophysical and geophysical systems. For instance, the A-effect 
and/or negative viscosity can provide a mechanism for the generation of ubiquitous large-scale shear flows, which 
are often observed in these objects. Furthermore, the anisotropic mixing of scalars should be taken into account in 
understanding the surface depletion of light elements in stars [45j |. Finally, we note that numerical confirmation of 
our prediction and the extension of our work to stratified rotating sheared turbulence with/without magnetic fields 
remain challenging important problems, and will be addressed in future publications. 
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APPENDIX A: SYMMETRIC PERTURBATION ((3 = 0) 

In this section, we consider a symmetric perturbation with k z = by assuming a forcing that is symmetric in the 
span-wise direction with no dependence on z. Note that even though k z = 0, v z and v y are closely linked through 
rotation Q,x. For j3 = k z /k y — 0, the homogeneous part of the first equation in (fUi| becomes: 

+ n 2 T« x = o. (Ai) 



Or 



r i 



d T ({l + T 2 }i x ) 
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Solutions of the homogeneous problem are thus Bessel functions. Using the method of variation of parameters, we 
can then express the general solution of the first equation to (|14[) as: 



2VT 



T 



TO 



(A2) 



Here again, fig = |f2|, \ — sign(f2); and L np are defined by: 

K P (t,T) = r„[n Vi + 1 2 ] J p [n Vi + t 2 ] - j n [n Vi + ^}Y P [fi Q ^i + ^} . (A3) 

The second equation of system (fT4"f can then be used to obtain the other components of the velocity in the form: 



TO 



^ X Loo(^,T)-M^^/^T^L 1 o(^ ) r) 



(A4) 



and a similar expression for v v (t). We can now use equations (|A2[) and (|A4|) to calculate turbulence amplitude f ^A 1[) 
and transport (i jA 21 and SjX3|. Note that equations (|A2|) and (|A4[) are exact solutions valid for all values of Sl - 



1. Turbulence amplitude 

From equations (|A2|1 and (|A4|) , we can easily obtain the turbulence amplitude as: 

2 _ T } ir 2 nl 



{IK 



4(2tt) 3 ^ 



J d 3 kF(k)(l + a 2 ) [Xi(k) + X 2 (k)] , (A5) 



(vl) = {f^yj I d 3 kF(k)(l + a 2 ) [X 3 (k) + X 4 (k)] 



Here, for simplicity, we considered only an isotropic forcing, given by equation (I12|) . and defined the following integrals: 



+°° e -mQ{-r)-Q{a)] 



Xi(k) = / T — 2 [L i(a,T)]^T, (A6) 



1 + T Z 
1+T 2 



e -2f[Q(r)-Q(a)] % 



X 2 (k) = / [L n (a, r)] z dr, 



/■+oo 

A 3 (k) = / e -^«W-«[L„„( a ,r)fdr, 

•/ a 

X 4 (k) = / e- 2 ^M-^>][L 10 (a,T)] 2 dT. 

Here, L„ p 's are given by equation (|A3|) . We now consider the strong shear limit: £ = vk 2 /A <C 1. As both Bessel 
functions becomes as (1 + t 2 ) -1 / 4 (up to a trigonometric functions) for large r, the first two integrals converge as 
£ — > 0. Thus, it is sufficient to put £ = in Xi and X 2 in equation (|A6|) to obtain the leading order behavior for 
In comparison, the integrand of A3 and A4 behaves as 1/r for r ^> 1, giving a contribution of order ln£ to 
leading order. 

We now examine the turbulence amplitude in the large rotation limit: Qq ^> 1. To do so, we use the asymptotic 
behavior of the integrals (1A6[) derived in appendix IB 1 bl Using equations (|B5[) and (|B7[) in equation(|A5p. we obtain 
the leading order contribution of the turbulent amplitude given by Eq. (|27[) in the main text. 

2. Turbulent transport of momentum 

We now calculate the turbulent viscosity vt defined by (v x v y ) = —vt9 x Uq = vtA. From equations (|A2|) and (IA4p . 
we can derive the Reynolds stress in the case of an isotropic forcing: 

(^v)^~i^^l d3k ma + a 2 )[X 5 (k)+X 6 (k)} , (A7) 
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where, 

*s(k) 

Here, L np 's are again given by equation (|A3|) . Note that the expression for the transport of angular momentum 
[equation (|A7j) ] is the same as that of (v 2 .) [equation ljA5|> ] except for the multiplicative factor of — r. This is simply 
because, for (3 = 0, the incompressibility condition imposes v y = —tv x . By using the asymptotic behavior of Bessel 
functions for large argument, we see that the two integrals X§ and Xq in equation (|A8I ) can be evaluated in the 
strong shear limit by just putting £ = 0. Consequently, the turbulent viscosity is of order A~ 2 for any value of f2. 

In the large rotation limit (CIq ^> 1), we can estimate the integrals (| A8|) and obtain the turbulent viscosity given 
by Eq. (J2HJ) in the main text. 



+ oo T e -2£[Q(r)-Q(o)] 



1 + T 2 

T e -2f[0(r)-Q(o)] 



1 + T 2 



[L i(i, r)] 2 dr , 



[Lu(a,r)] dr 



(A8) 



3. Particles transport 

The fluctuating concentration of particles can be obtained by integration of the fluctuating velocities (|A2[) and (|A4|) 
[see equation (JTU))]. Then, the diagonal part of turbulent diffusivity can be obtained as: 

^ = 4(C)M' / ^(l + « 2 W)[^i(k)+P 2 (k)] , (A9) 
^ = 4(C)3% / rf 3 Ml + « 2 W)^3(k)+P4(k)] . 
Here, we defined integrals Pi which all have the following form: 

/•-j-oo rr 

Pi(k)=/ dre- 2 « [(?(T) - Q(a)1 P(r) / F^tjdt, (A10) 



for « = 1 to 4. The functions p(r)'s are defined by: 

Ps = h 00 (a, r) , F 4 = -L 10 (a, r) . 

In the large rotation limit (f2 3> 1), the p's are oscillating functions. Thus, to evaluate integrals (|A10p in the 
strong shear limit (£ <C 1), we can not simply put £ = in equation (|A10|) as is explained in the appendix IB 31 A 
careful analysis (see appendix IB 3|) then gives us Eq. (p?5|) of the main text in the limits of strong shear (£ <C 1) and 
large rotation (CIq 3> !)• 



APPENDIX B: ASYMPTOTIC EXPANSION OF INTEGRALS 



In mil CI we took a large shear limit (£ <C 1) and obtain equation (|A5p for the turbulence intensity, equation (|A7|) 
for the transport of angular momentum, and equation (|A9[) for the transport of particles in terms of integrals involving 
Bessel functions of an argument depending on the rotation. We here derive asymptotic behavior of these integrals to 
simplify our results. 
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1. Non Oscillating integrands 

For non oscillating integrands, it is sufficient to put £ = in the integrals to find the large shear limit (the resulting 
integral converges as £ — > 0). Here, we provide asymptotic behavior of the following integrals for small or large Sl : 

I-+0O 1 

Xr(k) = J ___[ Loi (^ T )] 2 dTj (Bl) 



/•-t-oo i 

X2 ( k ) = / T— ^[Ln (a, r)] 2 dr, 
X 6 (k) - / — ^ [Ln (a, r)] 2 rfr. 



Here L„ p 's are given by equation (|A3 



a. Small rotation limit (T2o <§C 1 j 

To calculate Jfi and X5, one can use the asymptotic expansion of the Bessel functions and readily obtain: 

4 r°° dr 2 



Xr, 



^niJ a (i + r 2 ) 2 

rdr 2 



- - arctan(a)) - — — ^ , (B2) 
■ z 1 + a z J 



7T^7 o ( 7 + r 2)2 ^2 (1 + fl2) 



If we apply the same strategy to the calculations of Xi and Xq, the resulting expression would not be integrable so 
we have to calculate it otherwise: 



2 f°° J 2 {x) 8 



7r\/rTa 2 " 7o 3tt 2 VTTo? ' 

2 p rJ^Vl + r 1 ) , 2 /""Jfts) 

A fi ~ — , / — — ^ dr - ' 



7rf2 \/l + a 2 A 1 + r 2 7rf2 Vl + a 2 Jo ^ 7rSl \/l + a 2 

b. Large rotation limit (Qq S> 1) 

Using the Bessel asymptotic behavior for large argument, we obtain the following formula for the first integral: 

4 [+°° cos 2 [Oq{ VTT^ - VT+r^}] 

2 f +0 ° 1 2 . a 



7 r 2 n 2 v / r+^ 2 "7a (i + r 2 )3/ 2 7T 2 n 2 VTT^ x VTT. 

ans similarly for the other three integrals. Finally, we obtain the following asymptotic behavior for the four integrals 

Xl ~ X 2 - (1 - , a ) , (B5) 

•^5 ~ Xq 



■K 2 n 2 1 + a 2 
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2. Logarithmic divergence 

As noticed in $A 11 there is a logarithmic divergence arising in the calculation of X3 and X4. We here calculate this 
divergence in the case of a fast oscillation. Following [25[, we change the integration variable from r to y = 2£r 3 /3, 
replace the Bessel function by the expression valid for large argument (£ <C 1), and then obtain the following, to 
leading order in £: 



* 3 (ki) 



dy 



e~ v dy 



£a 3 



(3^)2/3(2^)1/3^1 + (|) 2/3 



(B6) 



{cos[^l + (|) 2/3 - jFoKa)] - sin[^l+(|) 2/3 - J] JoKo)]} 3 . 

We see that as £ tends to zero, the integrand in equation (|B6[) becomes proportional to 1/y, giving a contribution of 
the order ln£. 

In the large rotation limit (fio 3* 1)> we replace the Bessel functions by their asymptotic behavior to obtain: 



X, 



e v dy 



'y\ 2 / 3 ^ 



e~ v dy 



1 + (|) 2/3 



(B7) 



7T 2 ft 2 vT + a 2 J (a3 3 



ln£ 



2£jr 
3 



2/3 



tt^vT 



to leading order in £ <C 1. Following the same analysis, we find the same asymptotic behavior for Xa. 



3. Oscillating integrands 

The calculation of the transport of particles involves the computation of double integrals of the type: 

P= < tre- 3 *W- Q WF(T) F(t)dt, (B8) 



where the functions F contains an oscillating functions. We here deriv e the a symptotic behavior of this integral 
with F(t) — fit) cos[f2o</>(i)] and the phase given by 4>{t) = \fl + a 2 — yl + t 2 . The difficulty associated with the 
calculation of such integral is the presence of a point of stationary phase t — where the integral cannot be done with 
an integration by part. 

For a > 0, the point of stationary phase is never reached and then, the first integral can be approximated, for 
fl > 1, as: 



f T VI + t 2 fM 

I(t) = / F(t) dt ~ ~ n ' sin n ^(r) • (B9) 



Using this approximation, P can be computed with the following result: 

/(q) 2 (l + q 2 ) 
~ 4[£(l + a 2 ) 2 + 2 a 2 ] ' 



(BIO) 



Note that the result is the same as in the perpendicular case where the integral defining the transport of particles 
does not involve any stationary point. 

For a < 0, the behavior of the integral J(r) is affected by the stationary point in the vicinity of r = 0. We can 
however find an approximation as: 

.^Msin[Q o 0(r)] ^r<-^, 
J(r) ~ { h + ct if |r| < ^= , (Bll) 

2/o-^™ S m[O o 0(r)] ifr>^=. 
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FIG. 3: Graph of the function I(r) with our approximation ()B11|I . The parameters are a — —1 and fio = 100. 



Here, Iq = \J ir/2tt f(0) cos [Oq0(O) — tt/4] is the value given by the stationary point and c = /(0) cos [f2o</>(0)] is 
obtained by Taylor expanding / in the vicinity of r = 0. Figure [3] shows the numerical computation of the integral 
compared to the approximation (|B11[) and shows an excellent agreement. Using equation (|B11[) . we obtain P as: 

" /W 2 ( 1 + « 2 ) , 2/ 2. (B12) 



4[C(l + a 2 ) 2 + ^a 2 ] 

The first contribution comes from the integration by part (and as the result is odd in r, the contributions from 
— l/VT^o and l/vT^o cancel out). The second contribution (of order fi^ 1 ) comes from the stationary point. Both 
contributions have to be kept as the first one can be important if \uoa\ <C 1. 

For a = 0, the stationary point counts twice as less, so the approximation becomes: 

J cr if < r < 

/(T) ~U-*^HPo<Kt)] i^>^ ° (B13) 

In that case, the contribution from the stationary point cancels out as 1(0) = 0. Therefore, for a — 0, the only 
contribution comes from the end point of the integration and is the same as for a > [see equation (|BL0|) ] . 
Performing the same procedure when F(t) = f(t) sm[H,o(f>(t)], we obtain the following result: 

/(a) 2 (l + a 2 ) , /(a) 2 (l + a 2 ) 



4[£(1 + a 2 ) 2 + ft 2 a 2 ] [4^(1 + a 2 ) 2 + fi 2 a 2 ] 



7T 



^/(0) 2 sin 2 [fi o <K0) - tt/4] J , (B15) 
the second line being present only if a < (i.e. when the point of stationary phase is reached). 
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